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ABSTRACT 

We derive kinematic distances to the 86 6.7 GHz methanol masers discovered 
in the Arecibo Methanol Maser Galactic Plane Survey. The systemic velocities 
of the sources were derived from 13 C0 ( J = 2 — 1), CS ( J = 5 — 4), and NH3 ob- 
servations made with the ARO Submillimeter Telescope, the APEX telescope, 
and the Effelsberg 100 m telescope, respectively. Kinematic distance ambiguities 
were resolved using HI self-absorption with HI data from the VLA Galactic Plane 
Survey. We observe roughly three times as many sources at the far distance com- 
pared to the near distance. The vertical distribution of the sources has a scale 
height of ~ 30 pc, and is much lower than that of the Galactic thin disk. We 
use the distances derived in this work to determine the luminosity function of 
6.7 GHz maser emission. The luminosity function has a peak at approximately 
10~ 6 L Q . Assuming that this luminosity function applies, the methanol maser 
population in the Large Magellanic Cloud and M33 is at least 4 and 14 times 
smaller, respectively, than in our Galaxy. 

Subject headings: masers — Galaxy: kinematics and dynamics — Galaxy: struc- 
ture — radio lines: ISM 



1. Introduction 

Methanol masers at 6.7 GHz are unique compared to their OH and H2O counterparts 
in that they appear to be exclusively associated with early phases of massive star formation. 
They are hence extremely useful tools to identify and study very young massive star forming 
regions. To date, over 800 sources have been detect ed through vari ous targeted and blind 



surveys of the Galactic plane (e.g. iGreen et al.ll2009l : compilation of IXu et al.ll2009bl ). 
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In spite of the extensive surveys to date for 6.7 GHz methanol masers, their luminosity 
function remains largely unknown. Since 6.7 GHz methan ol masers are c losely associated 
with mainline OH masers in massive star forming regions, ICaswell et al.l (119951 ) suggested 
that the luminosity function of methanol masers would be similar to that of OH masers, 
but with a scaling factor to acco unt for the fact tha t meth anol masers are on average much 
more intense than OH masers. Ivan der Walt et al.l (119961 ) use a probabilistic approach to 
assign distances and estimated the luminosity function to have a power-law behavior with 
an index around -2. Both these studies used only the peak flux density rather than the 
integrated flux for determining the lumin osity functions. The most recent study of the issue 
was carri ed out bvlPestalozzi et al.l (|20070 w ho used a compilation of all sources detected prior 
to 2005 (IPestalozzi. Minier. fc Booth! 120051 ). Carrying out a statistical analysis of the maser 
population and modeling the spatial distribution of the masers, the luminosity function was 
modeled by those authors as a power-law with sharp cutoffs and having an index between 
-1.5 and -2. 

There are two problems facing studies of the luminosity function. First, most studies use 
a catalog of methanol masers compiled from various surveys, many of them targeted (towards 
IRAS sources or OH masers), and with different sensitivities. However, unbiased searches 
have all shown that targeted searches, especia lly towards IRAS sources, underestimate the 
number of sources by a fac tor of 2 or greater (IPandian fc Goldsmith! 120071 ; ISzymczak et al. 



20021 ; lEllingsen et al.lll996l ). For the best estimate of the luminosity function, one should 



employ a blind survey, since it is possible to analyze the limitations such as completeness 
reliably. Among the several blind surveys to date, by far the most sensitive one is the Arecib o 
Methanol Maser Galactic Plane Survey (AMGPS; IPandian. Goldsmith, fc Deshpandell2007l ). 
The AMGPS covered an area of 18.2 square degrees between Galactic longitudes of 35° and 
54°, and detected a total of 86 sources. This survey has a 95% probability of detection 
at a peak flux density of 0.27 Jy although sources as weak as 0.13 Jy have been detected. 
This makes the AMGPS catalog ideal for determining the luminosity function of 6.7 GHz 
methanol masers, especially at faint luminosities, although the relatively small size of the 
sample and area covered by the survey results in large statistical uncertainties. 

A more formidable problem is the determination of distances to the sources. Measuring 
distances is an old and challenging problem in observational astrophysics. While trigono- 
metric parallax is the most reliable method for determining distances, it cannot be applied 
readily to a large sample of sources. The usual technique for Galactic sources is to use a 
Galactic rotation curve to determine kinematic distances. While kinematic distances can 



have significant errors at times (e.g. Ku et al.l 120061 ). recent work using Very Large Base- 
line Interferometry (VLBI) has discovered systematic proper motions in young massive star 
forming regions, which in principle can be taken into account to improve kinematic distance 
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estimates (IReid et al.ll2009l ) 



A second challenge in the use of kinematic distances arises from an ambiguity between 
two distances (a "near" distance and a "far" distance) for sources located within the so- 
lar circle in the first and fourth Galactic quadrants. The most popular method to resolve 
the kinematic distance ambiguity is to measure an absorption spectrum towards the source 
(which has to have associated continuum radiation), and compare the velocit ies of the ab 



sorpt ion lines with that of the source and the tangent point (see Fig. 1 of iKolpak et al. 



20031 for an illustration of this technique). Distance ambiguities ha ve been successful - 
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HI absorption at 21 cm has been used by iPandian. Momjian. fc Goldsmith! (120081 ) to 
resolve the distance ambiguity towards 34 6.7 GHz methanol masers that have either directly 
associated 21 cm continuum, or belong to a cluster harboring a 21 cm continuum source. 
However, most methanol masers do not have any detectable radio continuum, presumably 
due to the young age of the exciting massive young stellar object. Hence, absorption line 
experiments can be used to resolve distance ambiguities towards only a small sample of 6.7 
GHz methanol masers. 



Burton. Liszt. &: Bakerl (119781 ) and iLiszt. Burton, fc Banial (Il98ll ) discovered that sev- 
eral HI self-absorption features in HI maps of the Galactic plane correlated with CO emis- 
sion features, and hence hypothesized that HI self- absorption could be used to determine 
distances to molecular clouds. HI self-absorption arises from cold HI in the foreground 
absorbing warmer radiation of the background at the a specific radial velocity within the 
velocity range covered by the background gas. Hence, only molecular clouds at the near 
kinematic distance can display HI self-absorption, since at the far distance ther e is no back- 



groun d emission at the radial velocity of the cloud. The theoretical work of iFlynn et al. 



( 120041 ) shows that molecular clouds have enough opacity in cold HI to exhibit self-absorption 
against stron g 21 cm backg r ounds (suc h as in the Ga l actic plane). HI self-absorption has 
been used by iJackson et al. J2002I ) and iBusfield et al.l (120061 ) to resolve the distance ambi- 
guity towards molecular clouds and massive young stellar candidates. 

In this paper, we describe work to resolve the distance ambiguity towards 6.7 GHz 
methanol masers detected in the AMGPS using HI self-absorption. The sources were ob- 
served in CO ( J = 2 — 1) and NH 3 to determine their systemic velocities and the line profiles 
of thermally excited molecular emission. The maser emission of most sources have multiple 
emission components, which at times are spread over as much as 20 km s _1 . While the central 
velocity of the maser emission is usually within a few km s _1 of the systemic velocity, it can 
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at times be offset by more than 10 km s _1 , and is hence not as reliable as velocities derived 
from thermal molecular emission. With distances determined from HI self-absorption, we 
will then look at the distribution of sources in the Galaxy and derive the methanol maser 
luminosity function. This in turn will be used to compare the methanol maser population 
in our Galaxy with that in nearby external galaxies. 



2. Observations and Data Reduction 

The CO ( J = 2 — 1) observations were made in December 2006 using the 10 m Submil- 
limeter Telescope^ (SMT) at Mt. Graham operated by the Arizona Radio Observatory. We 
used the 1.3 mm receiver using ALMA single sideband mixers. The mixers make both side- 
bands available simultaneously, the sideband rejection being 18-21 dB. We used the setup 
with the upper sideband including the 12 CO (J = 2 — 1) line and the lower sideband including 
the 13 CO ( J = 2 — 1) line. For a backend, we used the filterbanks (FFBs) which have 1 GHz 
bandwidth with 1024 channels yielding a velocity coverage of ~ 1400 km s _1 and a velocity 
resolution of 1.4 km s" 1 . The system temperature ranged from 175 K under good weather 
conditions to 600 K at low elevations in moderate weather. Since carbon monoxide is very 
prevalent in the Galaxy, we used absolute posi tion switching, choos ing the off-positions from 



the BU-FCRAO Galactic Ring Survey (GRS; IJackson et al.ll2006l ). The on-source integra- 



tion time was 2 minutes per source. Pointing and focus checks were done on Uranus at 
regular intervals. 



The data were reduced using the CLASS software (http://www.iram.fr/IRAMFR/GILDAS) 
and standard procedures. A first order polynomial was used for baselining the spectra. The 
data were then imported into IDL, and the systemic velocity was determined through a Gaus- 
sian fit of the 13 CO line using the procedure "xgaussfit"B We focused on the 13 CO spectra 
since this isotopologue avoids the problems with blended lines from the larger linewidths, 
and self-absorbed lines that are more common in 12 CO. Although the velocity resolution is 
1.4 km s _1 , most sources were detected at high enough signal to noise ratio to have their 
line velocities measured to an accuracy of a fraction of a km s _1 . In most cases, the CO 
line velocity that is within 5 km s _1 of the 6.7 GHz methanol emission was chosen as the 
systemic velocity of the source. 

A few sources displayed two 13 CO peaks within 5 km s _1 of the maser emission, resulting 



1 The SMT is operated by the Arizona Radio Observatory, Steward Observatory, University of Arizona, 
with partial support from the Mt. Cuba Astronomical Foundation. 

2 From [http: / / fuse.pha.jhu.edu/ analysis/fuse Jdl_tools.html| - the FUSE IDL tools website. 
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in an ambiguity in the determination of the systemic velocity of the maser. The CS (J = 5—4) 
line was observed in these sources using the APEX telescope^. Th e observations were c arried 
out in September 2008 using the APEX-1 single sideband receiver (jVassilev et al.ll2008l ). The 
receiver was tuned to the CS ( J = 5 — 4) line at 244.936 GHz. The Fast Fourier Transform 
Spectrometer (FFTS; iKlein et alJliooS ) was used for the backend with a bandwidth of 1 GHz 
and 16384 channels, yielding an effective velocity resolution of 0.12 km s _1 . The observations 
were made in position switched mode, with an on-source integration time of 3 minutes per 
source. The pointing was checked every hour, and the focus was verified at the start of 
the observations. These data were also reduced using the CLASS software. A low order 
polynomial baseline was subtracted and the center of emission was determined using the 
CLASS procedure "gauss" ("minimize" in CLASS90). 

We observed 15 sources for CS (J = 5 — 4) emission, of which 8 sources did not yield 
any detection to a la main beam temperature limit of ~ 0.04 K (at 0.6 km s -1 resolution). 
For these sources, we use the NH3 (1,1) transition observed at the 100 m MPIfR Effelsberg 
telescope^ (Pandian et al., in preparation) to discriminate between the two 13 CO velocities. 



3. Resolving the kinematic distance ambiguity 

The kinematic distance ambiguity was resolved by looking for the presence or absence of 
HI self-absorption at the systemic velocity of the source, with the presence of self-absorption 
implying the ne ar kinematic dis tance. The HI data were taken from the VLA Galactic Plane 



Survey (VGPS; iStil et al.ll2006l ). which has a resolution of V x 1' x 1.56km s _1 . Some sources 
are also close to compact HII regions which have had their distance ambiguity resolved using 
previous absorption line studies. In these cases, we verified the classification using HI self- 
absorption, and adopted the results from previous work in cases of ambiguous HI line profiles. 
Attention was paid to the line profiles of HI and 13 CO, which should be similar in cases of 
self-absorption. We also allowed for offsets up to ~ 1 km s~ x between the CO and HI lines 
when considering the presence or absence of self-absorption. The resulting distances and 
isotropic luminosities of the sources are shown in Table 1, and selected individual sources 
are discussed in section I3TT1 We show overlays of the HI and 13 CO spectra in Fig. [TJ 



3 This publication is based in part on data acquired with the Atacama Pathfinder Experiment (APEX). 
APEX is a collaboration between the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 

4 Based on observations with the 100-m telescope of the MPIfR (Max-Planck- Institut fur Radioastronomie) 
at Effelsberg. 



- 6- 



The distances in Table 1 are calculated using the modified kinematic distance technique 
outlined in Reid et al. (2009), which was found to minimize the error between trigonometric 
parallax distances and kinematic distances towards an ensemble of massive star forming 
regions. First, systemic velocities were calculated in a modified LSR frame of reference 
with (Uq, Vq, Wo) = (10.0, 5.25, 7.17) for the solar motion. Kinematic distances were then 
calculated using a flat rotation curve with Rq = 8.4 kpc, and Go = 254 km s -1 , with the 
assumption that massive star forming regions rotate 15 km s _1 slower than the circular 
rotation speed. Uncertainties were calculated by including a 7 km s _1 uncertainty in the 
modified LSR velocity. 



3.1. Notes on selected sources 

34-82+0.35- This source is close to the G34.4 massive star forming region (Galactic longitude 
and latitude coordinates ~ 34.4°, 0.23°). Its systemic velocity of 57.4 km s _1 , which is close 
to that of G34.4 (57 km s" 1 ), suggests that the two regions are physically associated. The 
HI spectrum at the source location does not show self-absorption, although nearby locations 
in the same cloud complex show prominent self-absorption. Since the G34.4 complex is 
thought to be at the near kinematic distance, this maser source is also tentatively placed at 
the near distance. 

35.03+0.35 - The HI line profile in t his so urce is ambiguous with an emission line inside 



an absorption feature. IWatson et al.l (120031 ) determine this source to be at the far distance 



using 6 cm formaldehyde absorption. Hence, this source is classified at the far distance. 

35.39+0.02 & 35.40+0.03 - These two sources are part of a cluster with the same CO 
velocity. There is weak HI self-absorption in both sources (only 35.40+0.03 is shown in 
Fig. 1) and hence both sources are classified to be at the near distance. 

35.59+0.06- There are two 13 CO lines in this source at 49.1 and 59.3 km s _1 . APEX CS 
observations show the 49.1 km s~ x component to be associated with the methanol maser. 
The lack of HI self-absorption at this velo city places thi s sourc e at the far distance. This 



is confirmed by a similar classification by iKolpak et al.l (120031 ) for a related source using 
HI absorption. 

35.79-0.17 - There are two 13 CO lines in this source at 57.5 and 62.4 km s _1 . There is 
CS emission at 62.5 km s _1 , while the NH3 line peaks at 61.3 km s _1 . Thus the 13 CO line 
profile arises from the presence of different cores in the molecular cloud rather than self- 
absorption. The presence of HI self-absorption near the CO velocity places this source at 
the near distance. 
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36.64-0.21 - There are two 13 CO lines at 75.2 and 80.8 km s -1 . The HI spectra show the 
75.2 km s _1 component to be at the far distance, while the 80.8 km s _1 component is at the 
near distance. The two features are also seen to be distinct clouds from the GRS datacube 
which shows the 80.8 km s _1 component to have a large angular extent, in contrast to 
the 75.2 km s _1 component which is extremely compact. The NH3 line shows emission at 
74.8 km s _1 indicating that the maser source is at the far distance. 

36.84-0.02 - The CO (J = 2 — 1) spectrum was taken at an incorrect location on the sky, 
and hence the GRS 13 CO (J = 1 — 0) spectrum is shown in Fig. 1 (the systemic velocity 
being verified by APEX CS observations). The presence of HI self- absorption places this 
source at the near distance. The GRS datacube also shows this source to be distinct from 
the source 36.70+0.09 even though they are at very similar velocities. 

37. 02-0. 03 & 37. 04-0. 04 - The source 37.04-0.04 ( J2000 coordinates: I8 h 59 m 4 2 .41, 3° 38' 32".8; 
peak flux density = 4.13 Jy) is a new detection discovered in MERLIN observations to de- 
termine accurate positions of the AMGPS sources (Pandian et al., in preparation). The 
GRS datacube shows both sources to be part of a molecular cloud that spans ~ 0.4° in 
Galactic longitude. The presence of HI self-absorption suggests that this cloud (and both 
maser sources) are at the near kinematic distance (only 37.02-0.03 is shown in Fig. 1). 

37.38-0.09 - There are two 13 CO lines in this source at 57.7 and 69.2 km s" 1 . NH3 observa- 
tions at Effelsberg show the 57.7 km s _1 component to be associated with the maser. The lack 
of HI self-absorption shows that this source is likely to be at the far dista nce. This agrees 



with t he distance determination of the nearby HII region 37.37-0.07 by iKuchar fc Bania 



(119941 ) based on HI absorption (towards continuum in the HII region). 



37.53-0.11 - The HI spectrum in this source shows several features with intensities below 
K. Interpreting these as absorption lines towards continuum emission, this source is at the 
far distance since absorption lines are seen between the systemic velocity of the source and 
the tangent point velocity (84 km s _1 ). 

37.55+0.19 - The 13 CO line profile suggests that this line is self-absorbed, which is con- 
firmed by Effelsberg NH 3 observations. Taking this into account, there is prominent HI self- 
absorption in this source, and hence we classify it at the near kinematic distance. 

37. 60+0. 42 - The HI spectrum shows absorption at the CO velocity, though the line profiles 
do not match well. The discrepancy however is not severe such as in 38.92-0.36 (below). 
Hence, we tentatively classify this source to be at the near distance. 

37 .77-0.22 - The HI spectrum in this source is similar to that in 37.53-0.11. The presence 
of absorption lines at velocities between the source velocity and that of the tangent point 
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indicates that this source is at the far distance. 

38.26-0.08 - There are two 13 CO lines at 11.6 and 17.7 km s _1 , the former being associated 
with the maser as deduced from NH 3 observations. The HI spectrum shows an absorption 
feature that is misaligned with the CO velocity by almost 2 km s _1 making it unlikely to be 
self-absorption. We thus place this source at the far distance. 

38.92-0.36- The HI spectrum shows a very broad absorption feature around the CO velocity. 
Since the feature is too broad to be interpreted as self-absorption from cold HI, and the CO 
and HI line profiles are different, we locate this source to be at the far distance, and attribute 
the HI feature to the intrinsic line profile rather than absorption. 

39.54-0.38 - This is another ambiguous case with the HI spectrum showing a broad absorp- 
tion feature that does not match the line profile of both 13 CO ( J = 2 — 1) and the ( J = 1 — 0) 
lines. Inspection of the GRS datacube suggests that this source is distinct from 39.39-0.14 
that is at similar LSR velocity (and which shows prominent HI self-absorption). Moreover, 
examining the HI and 13 CO (J = 1 — 0) spectra at other locations in the molecular cloud 
shows no self-absorption. Hence, this source is tentatively classified to be at the far distance. 

41-08-0.13 - The HI spectrum shows features that look like self-absorption at some loca- 
tions in the molecular cloud, while other locations including the maser position show no 
self-absorption. Thus, this is another ambiguous case whose distance ambiguity cannot be 
reliably resolved using HI self-absorption. Taking into account the results of the nearby 
sources at similar LSR velocities, 41.12-0.22 (the distance ambiguity for which has been re- 
solved by previous absorption line studies), 41.16-0.20 and 41.23-0.20, we place this source 
at the far distance though with low confidence. 

41.12-0.11 - The 13 CO spectrum in this source is complex with multiple emission compo- 
nents. Effelsberg NH3 data and APEX CS data show the 38.1 km s _1 component (which 
is also the strongest) to be associated with the maser. This HI spectrum shows absorption 
that is misaligned relative to the CO emission by more than 2 km s" 1 . Hence, we classify 
this source to be at the far distance. 



41.81 -0.10- Since the coordinates of this source are uncertain (as explained in lPandian. Goldsmith, fc Desh 



2007 1), we use the coordinates of the nearest 24 /im point source in the MIPSGAL survey 



( jCarey et al.ll2009i ). and also examine other locations in the parent molecular cloud. As in 
source 39.54-0.38, there is a broad HI absorption feature whose line profile does not corre- 
spond that of the relatively narrow CO line. Hence, this source is also classified to be at the 
far distance. 



42. 30-0. 30 - The HI spectrum shows an absorption feature that is offset relative to the CO 
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emission by ~ 2 km s 1 (similar to 38.26-0.08). Hence, this source is placed at the far 
distance. 

W49N region - There are 5 methanol masers in the W49N region: 43.15+0.02, 43.16+0.02, 
43.17+0.01, 43.17-0.00, and 43.18-0.01. The CO line shows emission at two velocities around 
and 10 km s _1 respectively. All sources in this region ar e placed at a distance of 11.4 ± 



1.2 kpc based on proper motion studies of H2O masers by iGwinn. Moran. fc Reidl (119921 ). 
This is also in very close agreement with the kinematic far distance (11.4 - 12.0 kpc) of the 
sources. 

43. 80-0. 13 - There is a HI absorption feature in the vicinity of CO emission although there 
is poor corresponde nce between the HI and CO line profiles. Absorption line studies of 



Watson et al.l (120031 ) find this source to be at the far distance. This highlights the importance 
of the similar line profiles in the HI and CO spectra for the HI feature to be described as 
self-absorption. 

44-64-0-52 - This HI spectrum in this source is similar to that in 38.92-0.36, showing a 
broad "absorption" feature and a narrow CO line, making it difficult to interpret the former 
as self-absorption. Hence, this source is classified to be at the far distance. 

W51 complex- There are 10 methanol masers in the AMGPS (48.89-0.17, 48.90-0.27, 48.99- 
0.30, 49.35+0.41, 49.47-0.37, 49.48-0.40, 49.49-0.37, 49.49-0.39, 49.60-0.25 and 49.62-0.36) 
that can be attributed to the W51 molecular cloud complex. Among these, 48.90-0.27 , 



48.99-0.30, and 49.35+0.41 belong to the 68 km s 1 component (ICarpenter fc Sanderslll998l ) 



while the remaining sources are in the W51 giant m olecular cloud. Using the trigonometric 



parallax distance of W51 IRS2 by Ku et al.l (j2009al ). these sources are placed at a distance 
of 5.11^4 kpc. This is also close to the kinematic tangent point distance (5.5 kpc) of the 
sources. 

52.92+0.41 - The LSR velocity of this source is greater than the tangent point velocity. 
Hence, this source is classified to be at the tangent point. 

53.14+0.07 & 53.62+0.04 ~ These two sources are part of a cloud complex that spans a 
range of ~ 1.2° in Galactic longitude. Both sources show prominent HI self-absorption and 
are hence at the near distance. 
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4. Discussion 

4.1. Comparing maser and systemic velocities 

Fig. [2] shows the histogram of the difference between (a) the peak velocity of maser 
emission (v p ) and the systemic velocity (v s ) and (b) the central velocity of maser emission 
(v c , the arithmetic average of the minimum and maximum velocities of maser emission; all 
quantities are shown in Table 1). The sources in W49N and W51 main are not included in this 
plot since the systemic velocities of these sources are not well known. Fig. [2] shows a puzzling 
double peaked histogram for the difference between the maser peak velocity and the systemic 
velocity. On the other hand, the central maser velocity has a much better agreement with the 
systemic velocity though with large scatter. The averag e value of v r — v. q is 0.08 ±3.8 km s~ l 



While this average value is similar to that obtained by ISzymczak. Bartkiewicz. fc Richards 



(120071 ). the scatter in the AMGPS sample is much larger. The fraction of sources where 
the difference betw een v c and v. q is less than 3 km s -1 i s also smaller in our sample - 68% 



versus 83% seen by lSzymczak. Bartkiewicz. fc Richards! (120071 ) in their work. Nevertheless 



Fig. [2] shows that in the absence of molecular line data, the central maser velocity is a better 
estimate of the systemic velocity than the peak emission velocity. 



4.2. Distribution in the Galaxy 

4-2.1. Galactocentric distance 

The top panel of Fig. [3] shows the distribution of Galactocentric distances for the 
AMGPS sample. It should be kept in mind that any survey that does not cover the en- 
tire Galaxy will have a selection function in terms of sampling the area/volume at different 
Galactocentric radii. The dotted line in the top panel of Fig. [3] shows an estimate of this 
selection function by calculating the fraction of the total area in each distance bin annulus 
that was within the survey limits (35° < I < 54°, where I is the Galactic longitude). A more 
unbiased description of this distribution is presented in the bottom panel of Fig. [3] where 
the surface density (in counts kpc -2 ) is shown as a function of Galactocentric distance. It 
should be noted that the survey limit of / = 35° translates to a minimum Galactocentric 
distance of 4.8 kpc, and so the surface density shown in the 4-5 kpc bin h as a large uncer 



tainty. The shape of the distribution is in very good agreement with that of iPestalozzi et al. 



(120071 ) although t he absolute numbers ar e higher by almost a factor of 3. Considering that 



the distribution of IPestalozzi et al.l (120071 ) is derived from a compilation of sources detected 



from a number of surveys (both targeted and unbiased) with different sensitivities, the close 
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agreement of the shape of the distribution is remarkable. The difference in absolute num- 
bers arises from both the shallower depth of previous surveys, and the fact that the entire 
Galaxy is not cover ed by the surveys in th e compilation (leading to the surface density being 
underestimated by iPestalozzi et al.l 120071 ). 



4-2.2. Locations in the Milky Way 

Fig. H] shows the locations of the AMGPS methanol masers in a face-on view of the 
Galaxy superposed on an artist's conception of the Milky Way. There is a significant cluster 
of sources near the Crux-Scutum arm (close to the Galactic bar). Taking the error bars 
into account, it is possible to reconcile the positions of most sources with those of the Crux- 
Scutum, Carina-Sagittarius and Perseus spiral arms. The sources 53.14+0.07 and 53.62+0.04 
are located in the Local arm, and 49.41+0.33 appears to be in the Outer arm. The three 
distant sources, 36.92+0.48, 38.66+0.08 and 42.70-0.15 lie in between the Outer arm and the 
Crux-Scutum arm, and may trace a spur from the Outer arm. Since most distances to the 
sources are kinematic and have significant uncertainties associated with them (even though 
an estimate of some of the p roper motions are taken into account by using the modified 



technique of lReid et al.ll2009i ). it is not possible to carry out any quantitative analysis of the 



spiral arms traced by AMGPS methanol masers. 



4-2.3. Vertical Distribution 



Fig. [5] shows the vertical distribution of the AMGPS sources about the Galactic plane. 
Also shown in Fig. [5] are Gaussian and exponential fits to the observed distribution. Both 
distributions peak at negative z, reflecting the fact that the Sun is located above the Galactic 
plane as defined by the IAU. The locations of the peak (13.1 pc and 12.8 pc for Gaussian 
and exponential fits respecti vely) are very close to th e estimated position of the Sun above 
the Galactic plane (~ 16 pc; lHammersley et al.lll995l ). The half- width at half maximum of 
the Gaussian fit is ~ 30 pc, while the scale height of the exponential fit is ~ 20 pc. These 
numbers are significantly less than the scale height of the Galactic thin disk, which is thought 
to be around 100 pc. This may in part be a selection effect caused by the relatively small 
coverage of AMGPS in ter ms of Galactic latitude (~ ±0.4°). However, a very similar scale 



height has been observed by Ivan der Walt et al.l (119961 ) using methanol masers detected prior 
to 1996. 



To investigate this issue further, we looked at the vertical distribution of infrared dark 



- 12 - 



clouds (IRDCs), which are thought to be sites of imminent and ongoing massive star for- 
mation. Since the clouds are seen as extinction features against the Galactic mid-infrared 
ba ckground, they are tho ught to be at t he ne ar kinematic distance. Using the catalogs 
of ISimon et al.l (120061 ) and iJackson et al.l (120081 ). which cover the first and fourth Galactic 
quadrants respectively, we determined the vertical distribution of IRDCs (Fig. [6]). The scale 
heights of Gaussian and exponential fits are essentially identical to what is seen for 6.7 GHz 
methanol masers. Moreove r, the scale heig ht of compact HII regions in the Galaxy is ob- 
served to be around 35 pc (IFish et al.ll2003l and references therein). This strongly suggests 
that the low scale heights observed for 6.7 GHz methanol masers using the AMGPS sample 
is real, and reflects the scale height of newly born massive stars in the Galaxy. 



4-2.4- Near/Far Asymmetry 



Of the 87 sources in the AMGPS, 21 are classified to be at the near distance, and one 
source is placed at the tangent point. We can estimate whether this is realistic by calculating 
the fraction of the volume sampled in the survey that is in between the Sun and the tangent 
point (and thus corresponds to the near kinematic distance). To get meaningful results, this 
must be weighted by some function that takes into account the lower density of sources at 
large Galactocentric distances. Using the surface density of Fig. [3] as the weighting function, 
28% of the sources shou ld be at the near dist ance. This is very close to the 24% observed in 
this study. In contrast, ISobolev et al.l (120051 ) suggest confinement of kinematic distances to 
the near kine matic distance to pro duce realistic statistical results. Our work coupled with 
the results of iKolpak et al.l (120031 ) . who o bserved ~ 25% of HI I regions to be at the near 



distance, indicates that the assumption of ISobolev et al.l (120051 ) that methanol masers are 
statistically likely to be at the near distance is flawed, and may be due to the fact that the 
masers considered by them were found in much shallower surveys than ours. Deep surveys, 
such as the AMGPS, are sensitive to most methanol masers in the Galaxy, which is reflected 
in the close agreement of the observed near/far asymmetry with that expected from volume 
considerations. 

A curious feature that can be observed in Table 1 is that all sources between Galactic 
longitudes of 40.3° and 53° are at the far distance. This may indicate a dearth of young 
massive star forming regions in the near side of the Carina-Sagittarius spiral arm between 
Galactic longitudes of 35° and 49° (corresponding to Galactocentric longitudes, (3, of 15° and 
45° in the terminology of an ideal log-periodic spiral arm, where f3 is measured clockwise 
from the line joining the Sun and Galactic center). Fewer than five sources in the AMGPS 
catalog have systemic velocities close to what is expected from this region (~ 25 km s" 1 at 
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I = 35° to ~ 55 k m s 1 at I = 49°). A similar paucity is seen in the IRDC catalog of 



Simon et al.l (120061 ). though there are a larger number of sources (in the radial velocity range 



given above) than in the AMGPS catalog. Accurate distance measurements to the candidate 
sources in the region, and future work with other star formation indicators is required to 
verify this observation. 



4.3. Luminosity function of 6.7 GHz methanol masers 

Once the kinematic distance ambiguity is resolved, one can calculate the isotropic lu- 
minosities of the masers (using their integrated fluxes) and the methanol maser luminosity 
function. The luminosities range from 4.0 x 10 -9 L to 2.0 x 10~ 4 L Q . The mean luminosity 
(computed logarithmically) is 8.5 x 10 -7 L & , while the median luminosity is 9.1 x 10~ 7 L & . 

Fig. [7] shows the luminosity function of 6.7 GHz methanol masers. Errors in the 
luminosity function may arise from three sources: (i) formal uncertainties in the calculated 
luminosities (which arise from uncertainties in the kinematic distance), (ii) the success rate 
of the HI self-absorption technique to resolve the kinematic dista nce ambiguity, and (iii) 



statistical uncertainties from the finite sample. iBusfield et al.l (120061 ) carried out a confidence 



test of the HI self-absorption technique by randomly shifting the CO spectrum in velocity and 
looking for self- absorption on the shifted data. Using this test, they estimate the success rate 
in determining the correct solution for the distance ambiguity to be ~ 80%. Hence, we carried 
out simulations in which 20% of the sources which had the distance ambiguity resolved solely 
from HI self- absorption had their resolution (near/far) reversed followed by recalculation 
of the luminosity function. The error bars indicated in Fig. [7| reflect the results of these 
simulations in addition to taking into account the formal uncertainties in the luminosities 
themselves, as well as the statistical uncertainties. 

To interpret Fig. [?[ we have to estimate the compl e tenes s of the survey for various 



luminosities. Fig. 3 of iPandian. Goldsmith. Sz Deshpande J2007h shows the completeness of 



the AMGPS as a function of peak flux density. However, the detection probability of the 
matched filtering algorithm used for source detect ion will be a function of the integrated flux . 



Hence, we repeated the simulations described in IPandian. Goldsmith, fc Deshpandd (120071 ) 
for different linewidths, and determined that the survey is complete for an integrated flux of 
0.08 Jy km s _1 . We can now calculate the limiting distance at which all sources at a given 
luminosity are detected. The fraction of the area defined by the limiting distance and the 
survey limits to the total area covered by the survey gives an estimate of the completeness at 
this luminosity. However, the surface density of methanol masers decreases with increasing 
Galactocentric radius as shown in Fig. [3J Hence, the area should be weighted by the surface 
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density to obtain an accurate estimate of the completeness. The dashed line in Fig. [7J shows 
the result of this calculation, indicating that all sources in the survey region with luminosities 
greater than 3.5 x 10~ 8 L & are detected by AMGPS. 

There are two caveats in the above discussion that should be noted. First, methanol 
masers are variable. However, significant variability (e.g. flaring behavior) is limited to a 
small number of sources. Since at any given point of time, one can expect some masers 
to be in a higher than normal state, while other masers would be in a lower than normal 
state, variability will not systematically affect the entire survey. The overall shape of the 
luminosity function would then be relatively unaffected when considering a large sample. 
Second, the luminosities calculated assume that maser emission is isotropic. On the other 
hand, maser theory shows that masers are beamed, with the nature and extent of b eaming 



being dependent on the geometry and whether the maser is saturated or not (e.g. lElitzur 



19921 ). However, little is known about maser g eometry or beaming in astro-physic al sources. 



Theoretical models of methanol masers (e.g. ICragg. Sobolev. fc Godfrey! 120051 ) assume a 
beaming factor, defined as the ratio of optical depths in radial and tangential directions 
(r esulting from an elongated g eome try along the line of sight), of 10. The VLBI observations 



of iMinier. Booth, fc Conway! (120021 ) resulted in the detection of core/halo structures, which 
could mostly be explained by saturated masers in a spherical geometry, although modeling 
implied large errors on the size of the halos. A spherical geometry however would not 
have any preferred direction for the beaming, and hence isotropic luminosities would still 
be accurate. Since the geometries of astrophysical masers are largely unknown, we cannot 
realistically address the effect of beaming on the luminosity function. 

Fig. [7] shows that the luminosity function rises sharply up to ~ 10~ 6 L & and declines at 
higher luminosities. Only the first two bins are severely affected by incompleteness (overall 
completeness ~ 95% for the third bin), and hence the continued decline in the number of 
sources for luminosities below 10 -7 L & appears to be real. At the high end of the luminosity 
function, we are limited by the relatively small coverage of the Galaxy by AMGPS. Since 
sources that have luminosities greater than 10 -4 L & are rare, one needs wide-area surveys 
to obtain an accurate census of such objects in the Milky Way. In any case, it is clear that 
the luminosity function cannot be described by a power law as was suggested in previous 
work. The behavior of the luminosity function for luminosities below ~ 10~ 8 L cannot be 
constrained by the AMGPS, or any other available sample. 

The luminosity function of 6.7 G Hz methanol maser s also appears to be different from 



that of OH masers as determined by ICaswell fc Hayned (119871 ). who found the number of 
OH masers to increase with decreasing luminosity down to 3 Jy kpc 2 (assuming that the 
Sun is located 10 kpc from the Galactic center). It should be noted that the unit Jy kpc 2 
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does not describe luminosity since the linewidth is not taken into account, and most OH 
masers, like Class-II methanol masers, display multiple emission components. Moreover, 
the kinematic distance ambiguity is not resolved, and the distinction between the two dis- 
tances is based on the expected number of sources at near and far distances from volume 
arguments. It is thus possible tha t the O H maser luminosity function is different from that 
determined by ICaswell fc Hayned (119871 ). However, constructing a figure similar to that of 
Caswell fc Havnesl Jl987h " using only the peak flux densities of the AMGPS sources results in 
a shape similar to that of Fig. [7] with a peak around 500 Jy kpc 2 with a decreasing number 
of sources at lower S p d 2 (where S p is the peak flux density and d is the distance). Hence, 
unless there are significant systematic errors in the distances to the sources, the luminosity 
functions of 6.7 GHz methanol masers and OH masers are dissimilar. 



The ongoing Methanol Multi-Beam survey (MMB; iGreen et al.l 120091 ) when completed 
should cover the entire Galactic Plane but at poore r sensitivity. The 90% completeness of 
the MMB is estimated to be between 0.9 and 1.1 Jy (IGreen et al.ll2009l ). which translates to 
an integrated flux of 0.34-0.41 Jy km s -1 (assuming a linewidth of 0.35 km s _1 ). A crude 
scaling from the sensitivity of the AMGPS would imply that the MMB will be complete to 
luminosities around (1.5 — 1.8) x 10~ 7 L , and will thus be able to verify the turn-over of 
the luminosity function. It is not clear whether the MMB will be able to probe the decline 
in luminosity function at lower luminosities, though it will determine the high-end very well. 



4.4. Extragalactic methanol masers 



We can use the luminosity function of methanol masers in our Galaxy to revisit the 
problem of dearth of such sources in external galaxies. To date, there h aye been searches for 
6.7 GHz methanol mas ers towards the Large Magellan i c Clo ud (LMC; ISinclair et al.lll992 



Ellingsen et al. 



1994 



(SMG lEllingsen et al 



Beas 



1994 



ev et al.l 



1996 



Beasley et al 



1996 



Green et al 



Green et al.l 120081) . th e Small Magellanic cloud 



2008h and M33 JGoldsmith. Pandian. fc Deslu 



20081 ). Four methanol masers have been discovered in the LMC, but there have been no de- 
tections in the SM C or M33. The 3<r sensitivity limit of the LMC and SMC surveys of 
Green et al.l (120081 ) are 0.27 Jy and 0.4 Jy respectively. In M33, the 3 sigma limit for an 
individual giant molecular cloud (GMC) is 4 mJ y, while the limit on the overall mean emis - 
sion by combining data for 14 GMCs is 1 mJy (IGoldsmith. Pandian. fc Deshpandd 120081 ). 
The distances adopted for LMC, SMC and M33 a re 50 kpc, 60 kpc and 730 kpc respectively 
dFeastlll999l : IWalkeri [l999l ; iBrunthaler et al.ll2005h . 



It must be borne in mind that while the luminosity function involves the total integrated 
flux, surveys are sensitive to the peak flux density of a source. This is especially a factor 
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for strong methanol masers which have a number of spectral features. Among the four 6.7 
GHz methanol masers in the LMC, two have a single spectral feature, while two have a more 



comp lex spectrum dominated by 3 features, two of which have similar strength (IGreen et al. 



20081 ). There are two approaches that one can use to evaluate the under- abundance of 6.7 
GHz methanol masers in external galaxies. The first is to use the peak flux densities of the 
AMGPS sources, scale them to the distance of each galaxy, and determine the number of 
sources that should have been detected in the surveys to date. A second approach would 
be to estimate a mean value of the integrated flux to peak flux density ratio, and use the 
luminosity function to estimate the number of sources above the threshold corresponding 
to the survey limits. In the AMGPS, the mean value of this ratio is close to unity (~ 1.3 
for very strong sources, and ~ 0.6 for weak sources). The results of the two approaches 
are tabulated under "Method-I" and "Method-IP respectively in Table 2. The number of 
methanol masers in the Galaxy has been assumed to be ~ 1300 (IPandian fc Goldsmithll2007l ; 
van der Waltll2005h . 



An additional factor to consider in this discussion is that the masses and star formation 
rates in the galaxies in question are different from that of the Milky Way, which will influenc e 
the total number of methanol masers in the galaxies. Using the estimates of [Israel (119801 ). 
the star formation rate in the LMC, SMC and M33 are 0.1, 0.02 and 0.2 times that of the 
Galaxy. If we assume that the number of young massive star forming regions is proportional 
to the current star formation rate (with an appropriate scaling for the total population of 6.7 
GHz methanol masers), then the total number of detectable methanol masers in the LMC, 
SMC and M33 reduce to 18-20, 2, and 14-18 respectively. The non-detection of methanol 
masers in the SMC is then not inconsistent with a methanol maser population with similar 
properties as that of the Galactic sources, though the much lower metallicity casts doubt on 
this assumption. The LMC population is about a factor of 4- 5 under-abundant c ompared to 
that in our Galaxy, which is consistent with the estimates of I Green et al.l (120081 ). The M33 
population is under-abundant by more than a factor of 14 although there is some uncertainty 
in the latter since only 14 GMCs were sampled in M33. 

This work thus confirms previous results that the methanol maser population in the 
LMC and M33 are different from that in our Galaxy. For a detailed discussion on the 
po ssible reasons for t he d earth of methanol masers in these galaxi es, we refer the reader 
to IGreen et al.l (120081 ) and iGoldsmith. Pandian. fc Deshpandd (120081 ) . However, the greater 



under-abundance in M33 in spite of its higher metallicity compared to the LMC is a mystery. 
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5. Conclusions 

We used 13 C0 (J = 2 — 1) observations and VGPS HI data to resolve the kinematic 
distance ambiguity towards 6.7 GHz methanol masers discovered in the AMGPS. The dis- 
tribution of the surface density of methanol masers as a function of Galactocentric distance 
agrees very well with existing estimates in the literature, although we observe the absolute 
numbers to be more than a factor of 3 higher. The vertical distribution of the sources has 
a scale height that is ~ 3-5 times lower than that of the Galactic thin disk, presumably re- 
flecting the smaller scale height of newly born massive stars. The resolution of the distance 
ambiguity allowed us to construct a reliable estimate of the luminosity function. Its shape 
does not agree with that of a power law, but has a peak around ~ 1CT 6 L Q followed by a 
decline towards lower luminosities. The luminosity function of 6.7 GHz methanol masers also 
appears to be different from that of mainline OH masers. Using the luminosity function, we 
derive estimates for the abundance of methanol masers in the LMC, SMC and M33 compared 
to the Milky Way. We find the under-abundance in M33 to be a factor of 3 higher than that 
in the LMC in spite of its higher metallicity. Finally, the distribution of sources between 
near and far distances closely follows the respective volumes sampled by the survey, thus 
indicating that the assumption of the near kinematic distance on a statistical basis should 
be avoided. 
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Table 1. Distances to 6.7 GHz methanol masers discovered in AMGPS. 
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60.0 


N 


3.9 + 0.4 


37.74 


(4.0 + 0.8) 


x 10" 


6 


36.02-0.20 


92.4, 


93.5 


93.0 


87.3 


N 


r 4 +0.8 


0.08 


i-6±8:i x 


io- 8 




36.64-0.21 


77.0, 


79.3 


77.3 


75.2 


F 


8.9 + 0.5 


0.52 


(2.9 + 0.4) 


x 10" 


7 


36.70+0.09 


52.2, 


63.2 


54.7 


59.7 


F 


9.7 + 0.4 


8.51 


(5.5 + 0.5) 


x 10- 


6 


36.84-0.02 


52.8, 


64.2 


61.7 


59.2 


N 


3.7 + 0.4 


4.65 


(4.4+1.0) 


x 10" 


7 


36.90-0.41 


83.1, 


85.1 


84.7 


79.7 


N 


5.0 + 0.6 


0.26 


(4.5 + 1.1) 


x 10" 


8 


36.92+0.48 


-36.3, 


-35.6 


-35.9 


-30.7 


F 


15.8 + 0.8 


0.51 


(8.8 + 0.9) 


x 10" 


7 


37.02-0.03 


77.5, 


85.3 


78.4 


80.6 


N 


5.1 + 0.6 


4.77 


8.611-j x 


io- 7 




37.04-0.04 


78.1, 


86.1 


84.7 


80.6 


N 


5.1 + 0.6 


6.66 


(1.2 + 0.3) 


x 10" 


6 


37.38-0.09 


67.5, 


70.9 


70.6 


57.7 


F 


9.7 + 0.4 


0.13 


(8.5 + 0.8) 


x 10" 


8 


37.47-0.11 


53.6, 


63.3 


54.7 


58.8 


F 


9.6 + 0.4 


22.41 


(1.4 + 0.2) 


x 10" 


5 


37.53-0.11 


48.2, 


56.6 


50.0 


53.0 


F 


9.9 + 0.4 


5.77 


(3.9 + 0.3) 


x 10" 


6 


37.55+0.19 


78.1, 


88.2 


83.7 


85.3 


N 


5-6t :6 


6.44 


1-418:! x 


10" 6 




37.60+0.42 


84.6, 


94.7 


85.8 


89.5 


N 


6-2±°;| 


25.29 


6.711:1 x 


10" 6 




37.74-0.12 


49.9, 


50.5 


50.3 


45.5 


F 


10.3 + 0.4 


0.25 


(1.8 + 0.2) 


x 10- 


7 


37.76-0.19 


54.9, 


66.0 


55.1 


59.2 


F 


9.5 + 0.4 


1.46 


(9.1 + 0.8) 


x 10" 


7 


37.77-0.22 


68.8, 


70.3 


69.6 


61.9 


F 


9.4 + 0.5 


0.37 


(2.3 ±0.3) 


x 10" 


7 


38.03-0.30 


54.6, 


65.9 


55.6 


62.0 


N 


3.9 + 0.5 


18.10 


(1.9 ±0.5) 


x 10" 


6 


38.08-0.27 


66.7, 


67.8 


67.5 


64.6 


N 


4.1 + 0.4 


0.18 


(2.1 ±0.4) 


x 10- 


8 


38.12-0.24 


66.5, 


79.7 


70.2 


83.2 


N 


°-°-0.6 


5.67 


1.218-j x 


10" 6 




38.20-0.08 


77.7, 


88.5 


79.6 


83.3 


F 


7 7 +0.6 
'■'-1.1 


16.26 


6.7lli x 


10" 6 




38.26-0.08 


6.1, 


15.9 


15.4 


11.6 


F 


12.2 + 0.5 


5.39 


(5.6 + 0.5) 


x 10" 


6 


38.26-0.20 


64.1, 


73.5 


70.2 


65.4 


F 


9.0 + 0.5 


1.28 


(7.2 ±0.8) 


x 10- 


7 


38.56+0.15 


23.1, 


31.2 


31.5 


29.3 


N 


2.1 + 0.5 


0.13 


4.0l?:° 7 x 


io-° 




38.60-0.21 


61.4, 


69.5 


62.6 


66.4 


N 


4.2 + 0.5 


0.59 


(7.2 ± 1.8) 


x 10" 


8 


38.66+0.08 


-31.9, 


-30.7 


-31.5 


-39.2 


F 


16.3 + 0.9 


0.71 


(1.3 ±0.2) 


x 10- 


6 


38.92-0.36 


30.8, 


33.5 


31.9 


38.7 


F 


10.5 + 0.4 


1.29 


(9.9 ±0.8) 


x 10- 


7 


39.39-0.14 


58.2, 


75.5 


60.4 


66.1 


N 


4.3 + 0.5 


0.77 


(9.9 ±2.3) 


x 10" 


8 


39.54-0.38 


47.4, 


49.4 


47.8 


60.9 


F 


9.0 + 0.5 


0.24 


(1.3 ±0.2) 


x 10" 


7 


40.28-0.22 


62.4, 


85.7 


73.9 


73.0 


N 


4 q+0-9 

^• a -0.6 


68.03 


1 1+ - 5 X 

1 - 1 -0.2 x 


10" 5 
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Table 1 — Continued 



Source 




^max ) 


v p 


v s 


KDA 


d 


s l 


L 






(km s- 1 ) 


(km s- 1 ) 


(km s -1 ) 




(kpc) 


Jy km s -1 


(L Q ) 




40.62-0.14 


29.7, 


36.7 


31.1 


32.3 


F 


10.5 + 0.4 


6.97 


(5.3 + 0.4) x 10- 


6 


40.94-0.04 


36.2, 


43.2 


36.6 


40.2 


F 


10.0 + 0.5 


1.50 


(1.0 + 0.1) x 10- 


6 


41.08-0.13 


57.2, 


58.4 


57.5 


63.8 


F 


8.4 + 0.6 


0.33 


(1.6 + 0.3) x 10" 


7 


41.12-0.11 


33.1, 


37.4 


36.6 


38.1 


F 


10.0 + 0.5 


0.63 


(4.4 + 0.5) x 10- 


7 


41.12-0.22 


55.0, 


66.6 


63.4 


60.1 


F 


8.7 + 0.6 


1.23 


(6.5 + 0.9) x 10- 


7 


41.16-0.20 


61.6, 


63.8 


63.6 


59.9 


F 


8.7 + 0.6 


0.20 


(1.0 + 0.2) x 10" 


7 


41.23-0.20 


54.0, 


64.9 


55.4 


59.2 


F 


8.7 + 0.5 


7.55 


(4.0 + 0.5) x 10" 


6 


41.27+0.37 


19.4, 


20.6 


20.3 


14.6 


F 


11.5 + 0.5 


0.16 


(1.5 + 0.2) x 10- 


7 


41.34-0.14 


6.6, 


15.0 


11.7 


12.6 


F 


11.6 + 0.5 


25.40 


(2.4 + 0.2) x 10- 


5 


41.58+0.04 


10.4, 


12.3 


11.9 


12.4 


F 


11.5 + 0.5 


0.22 


(2.0 + 0.2) x 10- 


7 


41.87-0.10 


15.5, 


23.7 


15.8 


18.8 


F 


11.1 + 0.5 


0.10 


(8.5 + 0.8) x 10- 


8 


42.03+0.19 


6.8, 


17.3 


12.8 


18.0 


F 


11.1 + 0.5 


30.10 


(2.6 + 0.3) x 10- 


5 


42.30-0.30 


26.2, 


34.7 


28.1 


27.1 


F 


10.5 + 0.5 


5.07 


(3.9 + 0.4) x 10" 


6 


42.43-0.26 


65.8, 


69.1 


66.8 


64.5 


F 


7.9 + 0.8 


1.61 


(7.0 + 1.4) x 10- 


7 


42.70-0.15 


-47.1, 


-39.0 


-42.9 


-44.3 


F 


15.9 + 0.8 


4.03 


(7.1 + 0.7) x 10- 


6 


43.04-0.46 


54.1, 


63.6 


54.8 


57.3 


F 


8.3 + 0.6 


10.39 


(5.0 + 0.7) x 10- 


6 


43.08-0.08 


9.6, 


14.9 


10.2 


12.3 


F 


11.2 + 0.5 


4.11 


(3.6 + 0.3) x 10- 


6 


43.15+0.02 


12.3, 


14.3 


13.3 


10.9 


F 


11.4+1.2 


12.03 


(1.1 + 0.2) x 10- 


5 


43.16+0.02 


6.8, 


22.1 


9.3 


3.5 


F 


11.4+1.2 


50.54 


(4.6 + 1.0) x 10- 


5 


43.17+0.01 


18.1, 


22.4 


19.0 


12.0 


F 


11.4+1.2 


15.05 


(1.4 + 0.3) x 10- 


5 


43.17-0.00 


-1.7 


, 4.2 


-1.2 


1.7 


F 


11.4+1.2 


1.44 


(1.3 + 0.3) x 10- 


6 


43.18-0.01 


10.3, 


11.6 


11.1 


11.9 


F 


11.4+1.2 


0.63 


(5.7+1.2) x 10- 


7 


43.80-0.13 


38.4, 


43.6 


39.6 


43.7 


F 


9.1 + 0.5 


43.18 


(2.5 + 0.3) x 10- 


5 


44.31+0.04 


55.0, 


56.6 


55.7 


57.0 


F 


7.9 + 0.7 


0.43 


(1.9 + 0.4) x 10- 


7 


44.64-0.52 


48.8, 


49.9 


49.3 


46.0 


F 


8.7 + 0.5 


0.29 


(1.5 + 0.2) x 10- 


7 


45.07+0.13 


56.8, 


60.0 


57.8 


59.2 


F 


7.4 + 0.8 


22.26 


(8.4+1.6) x 10- 


6 


45.44+0.07 


49.1, 


50.6 


50.0 


58.5 


F 


7.4±?; 6 5 


0.71 


2.7t?;4 x 10- 7 




45.47+0.05 


55.4, 


59.8 


56.0 


59.7 


F 




9.40 


3.4t?; 7 x 10" 6 




45.47+0.13 


57.0, 


73.5 


65.7 


61.9 


F 


6-9±S:g 


4.88 


(1.6 + 0.4) x 10- 


6 


45.49+0.13 


56.7, 


66.4 


57.2 


60.5 


F 


7 1+°- 7 

' ■ i -1.2 


4.98 


(1.7 + 0.5) x 10- 


6 


45.57-0.12 


1.2, 


9.8 


1.6 


4.8 


F 


11.2 + 0.5 


0.30 


(2.6 + 0.2) x 10- 


7 


45.81-0.36 


54.7, 


70.8 


59.9 


60.3 


F 


7.o±S:f 


8.41 


2.9±g;| x 10- 6 




46.07+0.22 


22.3, 


25.1 


23.3 


19.5 


F 


10.1 + 0.5 


1.20 


(8.5 + 0.8) x 10- 


7 


46.12+0.38 


57.5, 


62.9 


59.0 


55.0 


F 


' -°-i.o 


1.05 


4.lt?; 7 x 10- 7 




48.89-0.17 


57.2, 


57.5 


57.3 


55.7 


F 


5 1+ 2 -9 


0.02 


3.6 + f | x 10- 9 




48.90-0.27 


63.6, 


72.5 


72.0 


68.4 


F 


r i +2.9 
0J -1.4 


0.57 


i-o±S-J x io- 7 





Table 1 — Continued 



Source 




v p 


v s 


KDA 


d 


Si 


L 




(km s" 1 ) 


(km s _1 ) 


(km s- 1 ) 




(kpc) 


Jy km s -1 


(Lq) 


48. 99-0.30 


62.5, 72.6 


71.6 


67.5 


F 


5 1+2-9 
0.1_ 1A 


0.55 


l.OtH x 10- 7 


49.27+0.31 


-6.8, 7.5 


-3.2 


3.4 


F 


10.5 + 0.5 


13.71 


(1.0 + 0.1) x 10" 5 


49.35+0.41 


66.1, 69.2 


68.0 


65.5 


F 


5 1+2-9 
°- 1 -1.4 


6.72 


1.2+.J-J x 10- 6 


49.41+0.33 


27.0, -9.8 


12.1 


-21.3 


F 


12.2 ± 0.6 


23.18 


(2.4 ± 0.2) x 10 


49.47-0.37 


55.6, 76.1 


63.8 


62.4 


F 


5 1+ 2 -9 
0.1_ 1 4 


17.35 


3.1+tj x 10- 6 


49.48-0.40 


47.7, 65.3 


51.4 


58.2 


F 


5 1+2-9 
°- i -1.4 


14.90 


2.7tH x 


49.49-0.37 


54.4, 65.7 


56.1 


60.6 


F 


5 1+2-9 
°- 1 -1.4 


49.27 


0.9±H x 10" 5 


49.49-0.39 


49.9, 65.0 


59.3 


56.6 


F 


5 1+29 
°- i -1.4 


574.7 


1.01J-J x 10 ; 


49.60-0.25 


59.9, 66.7 


62.9 


57.2 


F 


5 1+2-9 


78.16 


1.4^-i x 10" 5 


49.62-0.36 


48.8, 60.1 


49.3 


54.9 


F 


5 1+29 


1.13 


2.0t^ x 10~ 7 


50.78+0.15 


47.6, 50.9 


49.1 


42.2 


F 


7.o±S:I 


2.79 


9.5±l;g x 10- 7 


52.92+0.41 


38.8, 45.0 


39.1 


44.8 


T 


5.1 


4.68 


8.4 x 10~ 7 


53.04+0.11 


9.7, 10.5 


10.1 


5.5 


F 


9.4 + 0.5 


0.69 


(4.2 + 0.5) x 10~ 7 


53.14+0.07 


23.4, 25.4 


24.6 


21.8 


N 


1.9 + 0.6 


0.54 


(1.4 + 0.9) x 10" 8 


53.62+0.04 


18.2, 19.5 


19.0 


23.3 


N 


2.1 + 0.6 


7.16 


2.2t\i x 10- 7 



Note. — The columns show the source name, range of maser emission (v m i n , v m ax), velocity of peak maser 
emission, v p , systemic velocity from molecular line measurements, v s , resolution of the KDA - near (N), far 
(F) or tangent point (T), distance, d, integrated flux, Si, and isotropic luminosity, L (calculated from Si). 
Asymmetric uncertainties in the luminosity are typically a consequence of the asymmetric uncertainties in the 
distance to the source. 



Table 2: Number of 6.7 GHz methanol masers that should have been detected in extragalactic 
surveys to date if the maser populations in the galaxies were similar to that in the Milky 
Way. 



Galaxy 


Limiting S p 
(mJy) 


Number of detectable masers 
Method-I Method-II 


LMC 


270 


180 


195 


SMC 


400 


100 


110 


M33 
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70 


90 
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Fig. 1. — Overlays of HI and 13 CO data. The solid line shows the HI spectrum derived 
from VGPS data at the maser position, while the dashed line shows the 13 CO spectrum at 
the same location. The 13 CO data correspond to the J = 2 — 1 transition unless otherwise 
mentioned in the text (§3.1). The left and right axes indicate the HI and 13 CO antenna 
temperatures respectively. The systemic velocity is indicated by the dotted line for clarity. 
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Fig. 1. — Continued. 
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Fig. 1. — Continued. 



-28- 




-20 20 40 60 80 100120 



1 20 
100 



20 






n Aaaa 40 - 62-0 - 1 4 
i V \ 


I 


i \ 

I \ - 







20 20 40 60 80 100120 

40.94-0.04 




0.8 

0.6 
0.4 
0.2 
0.0 



-20 20 40 60 80 100 120 





120 


-41.08-0.13 AA 




4 




100 






>> 


80 






3 


'U> 


60 






2 


c 

c 


40 


- / I 






c 


20 


- , .' 




1 















20 40 60 80 100 120 
-0.1 1 >\ 



20 20 40 60 80 100 120 

1 00 ^ 41.12-0.22 M ' « 

80 





-20 20 40 60 80 100 120 



-20 



20 40 60 80 100 120 

41.27 + 0.37^ 2 ' 5 




-20 20 40 60 80 100 120 

LSR veloc: 



-20 
ty (km/s) 



20 40 60 80 100 120 



Fig. 1. — Continued. 
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Fig. 1. — Continued. 



-30- 




-20 20 40 60 80 100120 



-20 20 40 60 80 100120 




-20 20 40 60 80 100 120 



45.47 + 0.13 



-20 20 40 60 80 100 120 





100 - 



/ i 


45.49 + 0.13 
I \ 1 


^ / ' 





20 40 60 80 100 120 



45.57-0.12 




-20 20 40 60 80 100 120 



45.81-0.36 




-20 20 40 60 80 100 120 

1 .5 



-20 20 40 60 80 100 120 



46.12 + 0.38 




-20 



20 40 60 80 100 120 -20 20 40 60 80 100 120 
LSR velocity (km/s) 



Fig. 1. — Continued. 
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Fig. 1. — Continued. 
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Fig. 2. — Histograms of the difference between the peak velocity of maser emission and 
the systemic velocity (top panel) and that of the mean velocity of maser emission and the 
systemic velocity (bottom panel). 
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Fig. 3. — The top panel shows the histogram of Galactocentric distances to the methanol 
masers in AMGPS. The dotted line (right ordinate) shows the fraction of area in the Galaxy 
of annuli represented by each distance bin that was covered by the survey. The bottom panel 
shows the derived surface density of methanol masers in the Galaxy. 



Fig. 4. — Locations of the methanol masers (blue dots with lines indicating the uncertainties) 
discovered in AMPGS. The Galactic center (red asterisk) and the Sun (red solar symbol) 
are located at (0, 0) and (0, 8.5) respectively. The background is an artist's conception of 
the Milky Way (R. Hurt: NASA/JPL-Caltech/SSC). The distances in Table 1 are scaled to 
an Rq of 8.5 kpc in this illustration. 
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Fig. 5. — Vertical distribution of the AMGPS sources. The dotted and dashed lines show 
Gaussian and exponential fits to the distribution respectively. The scale height derived from 
the fits are 30 pc and 20 pc from the Gaussian and exponential fits respectively. 
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Fig. 6. — The vertical distribution of Infrared Dark Clouds. As in Fig. O dotted and dashed 
lines show Gaussian and exponential fits to the distribution respectively. The scale heights 
derived for the infrared dark clouds are identical to that of 6.7 GHz methanol masers in Fig. 
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Fig. 7. — The luminosity function of 6.7 GHz methanol masers. The error bars shown include 
uncertainties in the kinematic distance, uncertainty in the resolution of kinematic distance 
ambiguity using HI self-absorption, and statistical uncertainties. The dotted line shows the 
completeness (right hand scale) of the AMGPS survey as a function of luminosity, the survey 
being essentially 100% complete for L > 3.3 x 10~ 8 L Q . 



